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We present a study of jet production in association with W and Z bosons in proton-proton collisions at a 
centre-of-mass energy of 7 TeV using the full 2010 data set collected by CMS corresponding to an integrated 
luminosity of (35.9 it 1.4) pb~-^. We report the measurement of ratios a(V + > n jets)/cr(V) and u(V + > 
n jets)/a-(V + > (n — 1) jets), where V represents either a W or a Z and n stands for number of jets. 



1. INTRODUCTION 

The study of the production of vector bosons W 
and Z in association with hadronic jets (known as "V 
+ jets" process), provides a stringent test of pertur- 
bative QCD calculations. Moreover, this process is a 
significant source of background in searches for new 
physics. Standard Model Higgs boson and for stud- 
ies of the top quark. A precise measurement of the 
cross section and an understanding of its kinematics 
is then essential. Finally, a clear signature associated 
at a relatively well known physics make this process 
an useful test for the commissioning of the detectors. 

We present the results of the analysis obtained 
with the full 2010 data sample collected by the Com- 
pact Muon Solenoid (CMS) experiment at the Large 
Hadron Collider (LHC) in proton-proton collisions at 
a centre-of-mass energy of 7 TeV and corresponding 
to an integrated luminosity of (35.9 ± 1.4) pb"^. To 
reduce systematic uncertainties, we measure the V -|- 
n jets cross sections relative to the inclusive W and 
Z cross sections. The complete V -I- jets analysis is 
reported in [T]. 



2. EVENT SELECTION 

A hardware-based trigger system selects a sample of 
electrons and muons then filtered in the online cluster 
with algorithms that have evolved following the rapid 
rise of the LHC luminosity. The largest sample of 
leptons preselected are electrons with pt > 17 GeV 
and muons with > 15 GeV. 

Monte Carlo (MC) simulation samples with a W 
or a Z boson are generated with the MadGraph event 
generator, producing parton-level events with a vector 
boson and up to four jets on the basis of a matrix- 
element calculation. MadGraph is interfaced to the 
PYTHIA program for parton shower simulation, modi- 
fied in such a way that the hardest emission is modeled 
using the exact matrix element calculation for one ad- 
ditional real emission. Top background processes are 
generated with MadGraph while QCD multi-jet and 



7 -I- jets processes are generated with PYTHIA alone. 
Whenever available, a NNLO or NLO normalization 
is performed. 

Minimum-bias events are superimposed to the hard 
interaction to simulate the event pile-up observed in 
data. 

A Particle Flow (PF) algorithm is used to recon- 
struct both the missing transverse energy {$t) ^^'^ 
the jets in the event. The PF algorithm creates a 
complete event description by collecting information 
from all of the sub-detectors and linking it together. 
Objects are formed into the categories of muons, elec- 
trons, photons, charged hadrons and neutral hadrons. 

Signal selection begins with the identification of a 
so-called "leading lepton" following the standard es- 
tablished by the measurement of the inclusive W and 
Z cross sections 0. 

For electron candidates we require a transverse mo- 
mentum selection px > 20 GeV and that the posi- 
tion of the ECAL supercluster lies in a fiducial region 
identified by |?7| < 2.5 with 1.4442 < \r]\ < 1.566 ex- 
cluded (this allow to reject electrons close to the bar- 
rel/endcap region transition where cables and services 
compromise the reconstruction). The selected elec- 
tron must match the object that triggered the event 
readout and it must also pass tight quality require- 
ments (including identification, isolation, and con- 
version rejection) which correspond to a lepton effi- 
ciency of roughly 80% as evaluated with a MadGraph 
-I- PYTHIA simulated sample. If a second electron of 
Pt > 10 GeV, detected within the ECAL fiducial re- 
gion, passes a looser set of quality cuts (corresponding 
to an efficiency of about 95%) and its invariant mass 
with the leading electron is in the range [60—120] GeV, 
then the event is assigned to the Z -t- jets sample. Oth- 
erwise, the event is assigned to the W -I- jets sample. 
Events with a muon with pT > 15 GeV and \ri\ < 2.4 
are rejected from the W -I- jets sample to reduce tt 
contamination. 

For muon candidates we require the presence of an 
isolated high quality muon with pT > 20 GeV in the 
region \ri\ < 2.1 with a transverse impact parame- 
ter \dxy\ < 2 mm to suppress cosmic-ray muon back- 
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Figure 1: Distributions of the Et for tiie leading jet (left) and of the number of reconstructed jets (right) in events 
W — >■ ^v. The ratio between the data and the simulation is also shown. The line at pr = 30 GeV corresponds to the 
threshold imposed for counting jets. 



ground. The isolation requirement is obtained by re- 
quiring that the combined activity of the tracker and 
calorimeters around the muon relative to the muon px 
is less than 0.15. If a second muon oi px > 10 GeV is 
detected in the region \ri\ < 2.5 such that the di-muon 
invariant mass lies in the range [60 — 120] GeV, then 
the event is assigned to the Z + jets sample. Other- 
wise, the event is assigned to the W + jets sample 

For the W + jets samples, the transverse mass 
Mt is computed starting from the lepton missing 
transverse energy using the relation Mt = 
^2pt-^t(1 ~ cos A0), where A0 is the angle on the 
plane orthogonal to the direction of the beams be- 
tween the direction of the pr lepton and the one. 
To avoid a region at low Mt containing essentially no 
signal we require that Mt > 20 GeV. 

Jets are reconstructed from the particle collection 
created with the PF algorithm and are formed with 
the anti-fcy clustering algorithm [5] with a size pa- 
rameter of R — 0.5. Jet energy corrections (JEC) are 
applied to flatten the jet energy response as a function 
of T] and pt- Jets must first satisfy identification cri- 
teria to eliminate those originating from noise in the 
calorimeter. We requested |?7| < 2.4 so that the jets 
fall within the tracker acceptance and, for the W sam- 
ple, a Mt > 50 GeV selection to reduce backgrounds. 

Selected events are assigned to exclusive bins of jet 
multiplicity by counting the number of jets with pT > 
30 GeV. The observed distributions for the leading jet 
transverse momentum and for the exclusive number 
of reconstructed jets in the W samples are shown in 
Figure [l] in data and simulation. The data is in good 
agreement with the MadGraph predictions normalized 
to the NNLO cross sections. 



3. SIGNAL EXTRACTION 

In order to provide model-independent results, we 
quote all results within the lepton and jet accep- 
tance, and only correct for efficiency of the selec- 
tion. The lepton efficiencies calculation is obtained 
from Z/7* -I- jets data samples by means of a tag- 
and-probe method, starting from selected events with 
two lepton candidates of invariant mass in the range 
[60 — 120] GeV. One lepton candidate, called the 
"tag", satisfies all selection requirements. The other 
one, called the "probe" , is selected with criteria that 
depend on the efficiency being measured. The signal 
yields are for two exclusive subsamples of events in 
which the probe lepton passes or in which it fails the 
selection criteria considered. Fits are performed to the 
invariant-mass distributions of the pass and fail sub- 
samples to extract the Z signal events. The measured 
efficiency is calculated from the relative level of signal 
in the pass and fail subsamples. The lepton selec- 
tion efficiency is the product of the reconstruction ef- 
ficiency, the identification and isolation efficiency, and 
the trigger efficiency, each of which is calculated as a 
function of the jet multiplicity in the event. 

For electrons, the efficiencies are roughly 70% (60%) 
for the W + jets (Z/7* + jets) signal events. 

For muons, the average efficiency is close to 85% 
(86%) for the W + jets (Z/7* + jets) signal events 
and it exhibits a dependence on the jet multiplicity. 

The signal yield is estimated using an extended like- 
lihood fit to the di-lepton invariant mass M;+;- for the 
Z + jets sample and to AIt for the W + jets sample. 

For the Z event samples, the main background pro- 
cesses, dominated by tt and W -I- jets, have little in- 
fluence and do not produce a peak in the Mi+i- dis- 
tribution. In the Mi+i- distribution we can then dis- 
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Figure 2: The ratio a(W + > n jets)/(T(W) and (t(W + > n jets)/(j(W + > (n — 1) jets) in the electron channel 
(left) and muon channel (right) compared to expectations from MadGraph (Z2 and D6T tune) and PYTHIA (Z2 tune). 



tinguish two components, one for the signal and one 
for all background processes. 

For the W sample, background contributions can 
be divided into two components, one which exhibits 
a peaking structure in Mt, dominated by tt, and an- 
other which does not, dominated by QCD multi-jet 
events. We perform a two-dimensional fit to the Mt 
distribution and the number of 6-jets, 'n!~^^^^'^^ . The 
Mt distribution allows the statistical separation of 
the signal from the non-peaking backgrounds, while 
^jot*'*^^°'^ distinguishes the signal and the other back- 
grounds from tt. 

In order to estimate the scaling rule of jets at the 
particle jet level, we apply an unfolding procedure that 
removes the effects of jet energy resolution and recon- 
struction efficiency. 

The main sources of systematic uncertainties in jet 
counting are the determination of the jet energy and 
jet energy resolution. We apply the jet energy scale 
(JES) corrections, available as a function of r] and 
Pt, to account for the detector response and inhomo- 
geneities. We apply also a pile-up energy correction to 
revise the pile-up subtraction method used that sys- 
tematically removes 500 MeV to jets in events without 
pile-up. While the systematic uncertainty in the jet 
counting is correlated among the different jet multi- 
plicities, the uncertainties from efficiency and fits are 
not. All statistical and systematic uncertainties are 
propagated in the unfolding procedure also estimat- 
ing the systematic error associated with it. 

4. RESULTS 

From the unfolded exclusive jet multiplicity distri- 
butions we derive inclusive jet multiplicities and cal- 



culate two sets of ratios: cr(V + > n jets)/(T(V) and 
a{Y + > n iets)/a(V + > (n- 1) jets), where CT(y) 
is the inclusive cross section. The results for the W 
sample are reported in Figure [2j where the systematic 
uncertainties associated with the JES and the unfold- 
ing are shown as error bands. For a large number 
of jets, the PYTHIA simulation alone fails to describe 
the data, while the MadGraph + PYTHIA simulation 
agrees well, as expected. 



5. SUMMARY 

We measured the rate of jet production in asso- 
ciation with a W or Z vector boson using the full 
2010 dataset collected by CMS in pp collision data 
at y/s = 7 TeV. We also measured the ratios of 
cross sections a(V + > n jets)/CT(V) and (t(V + > 
njets)/(T(V+ > (n — 1) jets) where n is the number 
of jets. All results agree well with simulations based 
on MadGraph -|- PYTHIA . 
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